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ABSTRACT 

Hie  flexibility  of  the  Fast  Field  Program  (FFP)  can  be 
improved  considerably  with  the  utilization  of  a technique  for 
collapsing  or  folding  the  kernel  of  the  function  to  be  Fourier 
transformed.  The  analysis  for  the  incorporation  of  this  method 
is  given  and  its  usefulness  is  discussed  and  demonstrated. 
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INTRODUCTION 

The  prediction  of  propagation  loss  at  close  ranges  poses  a 
severe  problem  for  normal  mode  calculations  at  both  ends  of  the  frequency 
spectrum.  This  predicament  is  related  to  the  number  of  discrete 
eigen  values  included  in  the  normal  mode  summation,  and  the 
exclusion  of  the  contribution  from  a possible  continuous  distribution 
of  eigen  values. 

The  range  to  which  the  contribution  from  the  continuous  spectrum 
is  significant  should  generally  decrease  with  increasing  frequency. 

The  number  of  discrete  eigen  values,  however,  generally  increases 
with  increasing  frequency  and  it  is  found  that  more  of  these  higher 
ordered  modes  are  needed  in  the  summation  to  achieve  accuracy 
at  close  ranges. 

As  the  frequency  decreases  the  concern  is  not  with  calculating 
enough  discrete  modes  but  with  the  possible  contribution  from  the 
continuous  portion  of  the  eigen  value  spectrum.  At  a sufficiently  low 
enough  frequency  it  is  possible  to  be  below  the  cut  off  frequency,  that 
is  no  discrete  modes  are  excited,  and  one  is  left  with  only  the 
continuous  spectrum.  Instances  of  this  type  will  be  the  subject  of 
a forthcoming  technical  note.  The  present  interest  is  with  the  case 
where  an  insufficient  number  of  discrete  modes  have  been  used  in 
the  normal  mode  summation. 

Although  the  Fast  Field  Program  (FFP)  is  not  an  eigen  value 
solution,  problems  can  arise  when  it  is  utilized  for  near  field 
calculations.  The  difficulty  lies  with  a sampled  wave  number  region 
which  does  not  include  a sufficient  portion  of  the  smaller  wave 
numbers.  This  situation  is  analogous  in  part  to  excluding  the  discrete 
higher  ordered  modes  in  the  normal  mode  summation.  FFP  calculations 
will  automatically  include  the  contribution  from  a continuous 
distribution  of  eigen  values  if  their  general  location  in  the  wave 
number  domain  is  included  in  the  sampling  region. 

Problems  of  a different  nature  occur  when  it  is  desired  to 
utilize  the  FFP  with  a so  called  "mini"  computer.  One  such  problem 
is  that  mini  computers  usually  have  a much  smaller  memory  than 
their  large  thi±d  generation  counterparts. 

In  this  note  a scheme,  uniquely  suited  to  the  FFP,  is  discussed 
which  provides  solutions  to  both  of  the  above  mentioned  difficulties. 
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THE  COLLAPSED  KERNEL 

The  pressure  field  for  a point  monochromatic  source  is  given 
(Eq.  (*0  of  reference  (l))  as 
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where 
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is  the  kernel  evaluated  at  the  discrete  values  of  the  wave  number’s 
horizontal  component ^ 
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The  field  is  given  at  the  discrete  ranges 
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provided  that  the  wave  number  sampling  distance,  d £ and  the 
range  resolution,  /j\.  are  related  according  to 
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where  N is  the  size  of  the  discrete  Fourier  transform. 

Before  discussing  the  Collapsed  Kernel  it  would  be  helpful 
to  emphasize  the  following  points: 

(a)  The  extent  of  the  sampled  wave  number  region  is 
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Following  the  procedure  given  by  A..  Hut tall  in  Section  2.5 
of  reference  (2)  it  ’.ill  be  shown  that  (l)  may  be  rewritten 
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where  the  Collapsed 


Kernel  ( -c  ) 
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In  comparing  (2)  to  (1)  we  find: 

(a)  The  extent  of  the  sampled  wave  number  region  is 
still  £ . 

To  see  this  consider  the  case  where  N=8  and  p=2.  Assume  that 
if  an  3 point  transform  were  done  as  in  (l)  the  magnitude  of  the 
kernel  is  given  as  shown  in  Sketch  1 
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Sketch  1 


The  Collapsed  Kernel  seen  by  Eq.  (3)  is  the  addition  of  the  magnitudes 
with  the  common  symbols  as  shovm  in  Sketch  2. 
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Sketch  2 


(b)  The  range  resolution  is  increased  to  — 

(c)  The  size  of  the  FIT  is  reduced  to  (A//f) 

The  integral  expression  for  the  field  is  given  approximately  by 
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jt  c'.ui  be  shown  that  beyond  same  point,/,  the  magnitude  of  the 
kernel  will  decay  exponentially.  Then  let  T be' given  by  which 

determines /V since  is  fixed  from  other  considerations. 

If  this  interval  of  integration  is  subdivided  into r segments 
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If  the  variables  are  evaluated  at  the  discrete  values 
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AN  EXAMPLE 

Propagation  loss  as  a function  of  range  at  200  Hz  is  shown 
in  Figure  1 for  a typical  summer  Mediterranean  profile.  With  the 
source  and  receiver  placed  near  the  surface,  50  and  100  ft  respectively, 
the  FFP  shows  a well  defined  image  interference  pattern  out  to  the 
first  convergence  zone.  In  obtaining  these  results  an  N=8l92  complex 
point  transform  was  utilized.  The  limits  of  the  sampled  wave  number 
region  were 

fc  --  . J 0 & o ?3  V©  = . JS3  3 7 £ 


and  the  distance  between  samples  was 
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The  results  are  correct  everywhere  except  in  the  range  interval 
from  about  two  to  ten  kyds.  This  error  can  be  traced  to  the 
fact  that  the  starting  wave  number,  , is  not  sufficiently 

small  enough.  The  problem  may  be  remedied  by  either  increasing  N, 
the  number  of  points  in  the  transform;  or  by  increasing  <3  £ , 
the  distance  between  samples.  The  former  would  require  a sustantial 
increase  in  computer  storage.  An  increase  in  the  sampling  distance 
could  lead  to  under sampling. 

A third  alternative  would  be  to  utilize  the  Collapsed  Kernel. 
The  Collapsed  FFP  prediction  for  this  example  is  shown  in  Figure  2. 

A ccmparison  of  the  two  predictions  reveals  that  collapsing  has 
altered  the  results  only  in  the  two  to  ten  kyd  range  interval  with 
the  collapsed  values  showing  the  correct  continuation  of  the  image 
interference  pattern.  This  was  acconplished  by  keeping  the 
same  and  A \ but  sampling  the  kernel  at  four  times 

as  many  points  so  that  the  new  starting  point  was 


The  32,767  samples  were  then  collapsed  four  times  according  to  (3) 
and  an  8192  complex  point  transform  was  performed.  In  this  instance 
collapsing  has  provided  the  correct  answers  by  increasing  the  length 
of  the  sampled  wave  number  interval  without  changing  the  memory 
requirements  or  the  sampling  distance. 
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Collapsing  could  also  be  of  value  ■when  it  is  desired  to  run  the 
FT?  on  a mini  computer.  Assume  that  for  some  case  it  is  determined 
that  an  8192  point  transform  is  needed  in  conjunction  with  some  <1  £ 
Also  assume  that  the  mini  can  only  handle  a 512  point  transform.  The 
FFP  could  be  run  on  the  mini  by  collapsing  the  kernel  sixteen  times. 
The  results  will  be  identical  to  those  obtained  by  doing  a straight 
8192  point  transform  except  that  <!/L  has  now  been  magnified  sixteen 
tines . 
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